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The idea to base the uncertainty relation for photons on the electromagnetic energy distribution 
in space enabled us to derive a sharp inequality that expresses the uncertainty relation [Phys. 
Rev. Lett. 108, 140401 (2012)]. An alternative version of the uncertainty relation derived in this 
paper is closer in spirit to the original Heisenberg relation because it employs the analog of the 
position operator for the photon — the center of the energy operator. The noncommutativity of the 
components of the center of the energy operator results in the increase of the bound 2>h/2 in the 
standard Heisenberg uncertainty relation in three dimensions. This difference diminishes with the 
increase of the photon energy. In the infinite-momentum frame, the lower bound in the Heisenberg 
uncertainty relations for photons is the same as in nonrelativistic quantum mechanics. A similar 
uncertainty relation is also derived for coherent photon beams. This relation has direct experimental 
consequences since it gives a precise relationship between the spectral composition of the laser beam 
and the minimal focal volume. 
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I. INTRODUCTION 

The nonexistence of the strictly localized photon states 
[1] and the associated lack of the photon position opera- 
tor makes it impossible to formulate the uncertainty rela- 
tion for photons in the standard Heisenberg form. On the 
other hand, it is obvious that also for photons the spread 
of momentum and the extension in space are subjected 
to some restrictions that embody the famous Heisenberg 
phrase [2] "Je genauer der Ort bestimmt ist, desto unge- 
nauer ist der Impuls bekannt und umgckchrt." 

Our approach to the photon uncertainty relations is 
based on two precisely defined concepts: the photon wave 
function in momentum space and the energy density of 
the quantized electromagnetic field. In our previous pub- 
lication [3] we used the second moment of the energy dis- 
tribution to measure the spread of the photon states in 
coordinate space. This led us to the uncertainty relation 
for photons in the form 



ArAp > Ah. 



(1) 



In the present work, we define the uncertainty of the 
position for photons that would be analogous to the stan- 
dard definition. For that we need some replacement for 
the (nonexistent) photon position operator. This role is 
played by R, which is the center of energy (or center 
of mass). The center of energy operator R is directly re- 
lated to the first moment of the energy distribution. This 
approach will allow us to obtain the uncertainty relation 
in a form even closer to the original Heisenberg relation, 



v / ar 2 Vap 2 > -n, 
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(2) 
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where d is the number of dimensions. A characteristic 
feature of the uncertainty relation for photons is that the 
left-hand side in this inequality in two and in three di- 
mensions is never equal to dH/2, but it tends to this limit 
with the increase of the average photon momentum. Only 
in the infinite-momentum frame is the uncertainty rela- 
tion for photons the same as for nonrelativistic massive 
particles. However, in one dimension, the inequality (2) 
is saturated so that in this case there is no difference 
between photons and massive nonrelativistic particles. 

We also prove the following sharp inequality: 



(3) 



In nonrelativistic quantum mechanics, the inequalities 
obeyed by the two measures of uncertainty, AR 2 AP 2 
and (R ■ R)(P ■ P), are completely equivalent. They 
have equal lower bounds and they are both saturated 
by Gaussian functions. This equivalence does not hold 
for photons. Nevertheless, the two inequalities are in- 
timately related. We shall first prove (3) and then use 
the information about the photon states that saturate 
this inequality to elucidate the intricate properties of the 
inequality (2). 

An early attempt to base the uncertainty relation 
for photons on the center of energy R was made by 
Schwinger [4], but he only gave a rough estimate that 
the lower bound of AR 2 AP 2 is of the order of ti 2 . 

In addition to an uncertainty relation for single pho- 
tons we derive a closely related uncertainty relation for 
photon beams. Using coherent states of the electromag- 
netic field to describe such beams in the limit of a large 
number of photons we prove the following sharp inequal- 
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ity: 
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VaWVap 2 >-h\ i + 
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(4) 



and we find the mode functions of the coherent states 
that saturate this inequality. 

II. THE CENTER OF ENERGY 



positive square root is unique. The significance of N is 
further underscored by its being the generator of Lorentz 
transformations. Since the operators H and N do not 
commute (the energy changes under Lorentz transforma- 
tions), the equivalence of the two forms of R in (9) is not 
obvious and is proved in Appendix A. 

It follows from the commutation relations between the 
generators of the Poincare group [7, 8] , 



[Ni 



ihSijH, 



(11) 



The nonexistence of the local photon density in config- 
uration space is due to the fact that in quantum electro- 
dynamics the operator of the total number of photons N 
involves not a single but a double integral [5]: 



N = 



1 



d 3 i 



d 3 r' 
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B(r,t)-B(r',t) 
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F^r,t)-F(r',t) 
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(5) 



We use systematically the Riemann-Silbcrstein vector 
(the RS vector) [6] 



~, . D(r,t) B(r,t) 



2fi 



(6) 



which will allow us to write many formulas in a com- 
pact form. The normal ordering removes the (infinite) 
contribution from the vacuum state. In contrast to the 
total-number operator, the total-energy operator of the 
electromagnetic field H (the Hamiltonian) is an integral 
of a local density, 



H= d 3 r£(r,t), 



where 



£(r,t) =:F^(r,t)-F(r,t): 



(7) 



(8) 



The center of the energy operator can be introduced in 
any relativistic theory. All we need for this construction 
is the set of generators of the Poincare group. Following 
Born and Infcld [7], we define the operator R as follows: 



R= -Ln + N-L 

2H 2H 



1 N 1 



H 



H 



(9) 



where N is the first moment of the energy distribution, 

N = jd 3 rrE{r,t). (10) 

The symmetrization in (9) is necessary to obtain a Hcrmi- 
tian operator. The inverse of the Hamiltonian is well de- 
fined, provided we exclude the vacuum state. The spec- 
trum of the Hamiltonian is nonnegative, therefore the 



that R and the total momentum P obey the canonical 
commutation relations between the position and momen- 
tum, 



[Ri,Pj] = ihSij. 



(12) 



We must, however, resist the temptation to treat R as 
a bona fide position operator because its components do 
not commute, 



[Ri, Rj] = -itupH^SijH- 1 , 



(13) 



where Sij is the operator of the intrinsic angular momen- 
tum: the difference between the total angular momentum 
and the orbital angular momentum, 



Si 



Mtj- 



R l P 1 - RjPi 



(14) 



Note that the effects of the noncommutativity are present 
in all systems with intrinsic angular momentum and de- 
crease with the increasing energy. We shall fully confirm 
this observation in Sec. VI. 



III. RELATIVISTIC UNCERTAINTY 
RELATIONS IN ONE, TWO, AND THREE 
DIMENSIONS 

Despite all of the differences between the nonrcla- 
tivistic and relativistic dynamics we may derive a sharp 
Heisenberg uncertainty relation along one direction, say 
x, for any relativistic system. This onc-dimcnsional un- 
certainty relation is based solely on the commutation re- 
lations between X = R x and P = P x and has the stan- 
dard form 



where 



AX 2 = ((AP) 2 ), 
AP 2 = ((AP) 2 ), 



AX 

AP: 



= X 
P- 



-(X), 

(P)- 



(15) 

(16a) 
(16b) 



The one-dimensional uncertainty relation holds for any 
relativistic quantum system. A simple proof of (15) uses 
the commutation relations (12) and the non-negative ex- 
pectation value of the operator: 



(AA-UAP) (Al + UAP)) >0, (17) 



)> 
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where A is an arbitrary real number. The condition that 
this expression treated as a function of A can have at 
most one real root gives (15). This inequality is satu- 
rated by the quantum state whose state vector satisfies 
the condition 



AX + iXAP) |*) = 



(18) 



The specific form of depends, of course, on the system 
under study. Note that we may remove the average values 
(X) and (P) from (18) by choosing |\l>) in the form 

|*) = cxp (i{P)X/h-i{X)P/h) |*')- (19) 

Since the inequality must hold for all vectors, replacing 
by \SSf') makes no difference and the two forms of the 
uncertainty relation in one dimension, namely, 



Vax^VaF 1 > \h 



and ^¥)^)>\h, (20) 



are completely equivalent. In nonrelativistic quantum 
mechanics the equivalence holds in any number of di- 
mensions. A spherically symmetric Gaussian function 
shifted in the coordinate space by (r) and in the momen- 
tum space by (p) by the unitary transformation of the 
form (19) will automatically saturate the inequality (2). 
This equivalence, however, is no longer valid for relativis- 
tic systems in three dimensions. 

To extend our analysis to two and three dimensions, 
we introduce the dispersion in position that involves two 
or three components of the center-of-energy vector R, 



AR 2 = (AP-AP), 



(21) 



where AR — R — (R) and the dispersion in momentum, 

AP 2 = (AP-AP), (22) 

where AP = P — (P) . Following the same procedure as 
the one used in deriving (15), we obtain (2). The proof is 
based this time on the expectation value of the following 
positive operator: 

(AP - UAP) • (ar + iAAp) J> > 0. (23) 

In contrast to the one-dimensional case, the inequali- 
ties (2) and (23) are not sharp because there is no state 
vector that is annihilated by all three components of the 
vector operator A = AR + i\AP and even by two com- 
ponents. This is due to the fact that the commutators 
(13) of the components of R do not vanish. Should there 
exist a state vector annihilated by A, then this vector 
would also be annihilated by the commutators of the 
components of A. These commutators are proportional 
to the components of spin. Therefore, for any relativis- 
tic quantum system endowed with spin the inequality (2) 
cannot be saturated. 

In the next section, we introduce a convenient formal- 
ism to describe photon states that will be later applied to 
derive the inequalities (3) and (4) and also to elucidate 
the meaning of the inequality (2). 



IV. QUANTUM MECHANICS OF PHOTONS 

In what follows, we shall consider one-photon states 
of the electromagnetic field. These states are generated 
from the vacuum state by the action of the photon cre- 
ation operators, 

— [Mk)al(k) + .U(k)aUk)\ |0), (24) 

where a± (fc) create photons with momentum hk and pos- 
itive or negative helicity A (left-handed or right-handed 
circular polarization). We assume the normalization of 
these operators such that the commutation relations have 
the form 



a x (k),a[,(k')] =S xy kS (3 \k-k r ). 



(25) 



This leads to the relativistic form (the volume element 
on the light cone d 3 k/k is invariant under Lorentz trans- 
formations) of the scalar product, 

(/ (i) i/ (2) )=^^/r(fe)/f(fe), (26) 

and the associated norm of one-photon state vectors, 

</l/>HI/l| 2 =^l/A(fc)| 2 - (27) 

The symbol ^ stands for the summation over A and the 
integration over k, 



(28) 



The functions /+ (fc) and /_ (fc) are the photon wave func- 
tions in momentum space. Their moduli squared are the 
probability densities to find the left- or right-handed pho- 
tons with momentum hk. 

The creation and annihilation operators are connected 
with the field operators through the expansion of the RS 
operator into plane waves [9], 



F(r,t) = Vfc J - 



d 3 k 



(2tt) 3 / 2 

xe(fc) a+{ky k - r - lut + aUk)e- lk - r+tult ] . (29) 



The normalized polarization vector e(fc) is: 
fc x (n x fc) — ik(n x fc) 



e(fc) = 



y/2k\n x fc| 



(30) 



where n is an arbitrary unit vector. 

In order to find the action of all relevant operators 
on one-photon states, we first express these operators in 
terms of creation and annihilation operators. This task 
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is simplified by using the RS vector in the form (29) and 
we obtain [9] 



H = ^lfuva[(k)a x (k), 



— hk a\(k)a x (k), 



(31a) 
(31b) 



M = ^lhal(k)(kx- i D x + X^ja x (k), (31c) 

N = ^^huja{(k)iD x a x (k), (31d) 

where D\ is the covariant derivative in momentum space 
on the light cone, 



D\ = V — iXa(k), 



(32) 



a( fe) =^( fe) .V e(fe) = ^^, (33) 

the dot denotes the scalar product of polarization vectors, 
and V denotes the derivatives with respect to k. 

In relativistic quantum mechanics of photons, the gen- 
erators of the Poincare group (31) act on the photon wave 
functions as follows: 



Hf x (k) = fkjf x (k), 
Pf x (k) = hkf x (k), 

Mf x (k) = h(kx- i D x + X^)J x (k). 

Nf x (k) = hujiD x f x (k), 



(34a) 
(34b) 

(34c) 
(34d) 



where we stretched our notation keeping the same sym- 
bols to denote the operators acting on the states of the 
field and the operators acting on the photon wave func- 
tions. Since all of these operators are Hermitian with 
respect to the scalar product (26), they generate two 
unitary representations f+(k) and f-(k) of the Poincare 
group. These representations are concrete realizations of 
the general scheme described in [10]. 

The center-of-energy operator R given by the sec- 
ond expression in (9) has the following representation 
in quantum mechanics of photons: 



1 



Rfx{k)=iVkD x -j=f x ik). 



(35) 



It is often convenient to replace the function f x (k) by its 
rescaled counterpart g x (k), 



9x(k) 



fx(k) 



(36) 



The transformation properties of g x (fc) under the Lorentz 
transformations are more complicated than those of 
f x (k), but this function is similar to the nonrelativis- 
tic wave function because in contrast to (27) its norm 



(and also the scalar product) has a familiar nonrelativis- 
tic form 

N| 2 = ^>!sa. (37) 
The center-of-energy operator acting on g x (fc) is 



Rg x {k) = iD x9x (k) 



(38) 



As a simple application of this formula, we find now the 
function that saturates the general one-dimensional un- 
certainty relation (15) in the case of photons. Choosing 
the direction in this relation along the n vector, we find 
that the covariant derivative (32) becomes an ordinary 
derivative along this direction because the component of 
ot(k) along n vanishes. Therefore, the function g x (k) 
which saturates the inequality is a Gaussian in the direc- 
tion n. This result has been obtained before by Holcvo 
[11] in the framework of estimation theory. 

The extension of the Heisenberg uncertainty relation 
for photons from one to three dimensions is nontrivial. 
In the next section, we use the representation (38) of the 
operator R to fulfill this aim. 



V. UNCERTAINTY RELATION FOR THE 
PRODUCT OF (R R) AND (P P) 

The formulation of the uncertainty relation for the pho- 
ton will be carried out with the use of the operators R 
and P acting on the photon wave functions /A(fc) in mo- 
mentum space. In this section we shall consider the prod- 
uct of the quantities (R ■ R) and (P ■ P), instead of their 
variances. The variances AR 2 and AP 2 reduce to (R-R) 
and (P ■ P) only when both (R) and (P) vanish. 

The quantities (R ■ R) and (P ■ P) expressed in terms 
of the rescaled wave function g x are: 

^■ ft ) = ^f^ D ^xT-D xgx 

= ^^{VgyVg x + \ 2 a 2 (k)g x g x 



+ iXa(k)-(g* x Vg x -g x Vg* x ) 



(P-P) 



3a- 



(39) 



(40) 



There is one immediate conclusion that can be drawn 
by inspecting the integrand in the formula for (R ■ R) . 
Namely, the presence of ot(k) rules out spherically sym- 
metric functions. To obtain a finite value of (R ■ R) we 
must eliminate the singularity at \n x fc| = by the 
appropriate angular dependence of g x . Our analytic so- 
lution will confirm this expectation. The breaking of the 
spherical symmetry is an important difference between 



■5 



the uncertainty relation for photons and for the nonrela- 
tivistic particles. 

Further calculations are most easily done after the 
transformation of the integrals to spherical coordinates, 

(R-R) = t— V/ dkk 2 / d6 sine d<p 
\\9\\ x Jo Jo Jo 

\ dkgxl 2 + \de 9 x\ 2 , \d v9x \ 2 , A 2 cos 2 % A | 2 



+ 



(P.P) = 



k 2 k 2 sin 2 9 
i\cos9(g* x d v g x - g x d v g* x ) 



k 2 sin 9 



k 2 sin 2 9 



(41) 



fe2 rco r-T rl-K 

tt-thV dkk 2 dO sin(9 dtpk 2 \g x \ 2 , (42) 
\\9\\ ^ Jo Jo Jo 



\g\\ 2 = V/ dkk 2 \d9 sin 9 d^\g x \ 2 . (43) 
, Jo Jo Jo 



The left-hand side of the uncertainty relation for R ■ R 
and P ■ P divided by H 2 is a dimensionless quantity which 
will be denoted [12] by j 2 , 



7 



2 _ (R-R)(P-P) 

h 2 



(44) 



We shall determine the minimal value of 7 applying 
a variational procedure, as we have done in [3]. The 
variation of j 2 with respect to g x (k) leads to the following 
equation for g x (k): 



1 1 \ 2 

d K r 2 d K fa-ci^aa. 1 ^2 \2 ot\„™aa \ 

k 2 k 2 sin 2 1 



(d e sin 9dg + d 2 — X 2 — 2i\ C os9d v )- — + jk 2 - 2 7 



H 



gx(K,9,<p) = 0, 



(45) 



where we replaced k by the dimensionless variable n, 

1/4 



(P-P), 



(46) 



After performing the variation, we put ||g|| 2 = 1. The 
variational equations for two values of A decouple, so that 
we may take one value of A at a time. Since the change 
of the sign of A is compensated by complex conjugation, 
we will consider only positive helicity A = 1. 

Equation (45) allows for the separation of variables, 



g x {K,6,<p)=K{K)Q(6)e im *, 



(47) 



and we obtain the following equations for the radial and 
the angular parts: 



— ^d K K d K + 



i(i + i)-A 2 



+ 7« 



JC(k) - 2 7 /C(k), 
(48) 



sin0 



dg sin 9 do + 



m 2 + A 2 - 2Am cos ( 



e(9) 



sin 2 9 

= 3(3 + 1)6(0). (49) 



The equation for 9(6*) is the same as in the theory of 
magnetic monopolcs (cf. [3], [13]). Its solutions are given 

in terms of Jacobi polynomials pj m,m )(#) (also known 
as "monopole harmonics"), 

e(9) = sin A #cot m 0P^ m m ' A+m) (cos#). (50) 



Regular solutions are obtained when j is a natural num- 
ber starting from j = 1. For j = 0, not only are both 
solutions of the angular equation, namely, 1/ sin 9 and 
cot 9, singular, but also the radial equation does not have 
regular solutions because the centrifugal force becomes 
attractive. Therefore, the s states are ruled out as we 
already observed before. 

The equation for JC(k) is the radial part of the 
Schrodinger equation for the three-dimensional harmonic 
oscillator with a modified centrifugal force. This equa- 
tion after the substitution, 



£(«) 



-3/2 



exp (- 



¥ 2 ) 



(51) 



reduces to the equation for the confluent hypergeometric 
function, 



7 2 



(52) 



where v — 1 + yj j + j 2 — 3/4. To obtain a regular solu- 
tion, \F\ must become a polynomial and this leads to 
the quantization condition for the parameter 7, 



7 = 2n + 1 + y/j + j 2 - 3/4, n = 0,l,2.... (53) 
The lowest value of 7 is obtained for j = 1 and n = 0, 



7 




(54) 



This is the right-hand side in the uncertainty relation (3). 

In what follows, we shall denote by 7 always its low- 
est value (54). This eigenvalue is degenerate. There are 



6 



three eigenfunctions that saturate the inequality (3) cor- 
responding to to = 0, ±1: 

/o(fc,M) = Aasin6»(afc) 7 - 1 cxp(-i(afc) 2 ) , (55a) 

f ± (k, 9, <f>) = Aa ^^^ e^iaky- 1 exp (-±(afc) 2 ) , 

(55b) 



where the normalization constant is 



47rr( 7 ) ' 



(56) 



and the parameter a sets the length scale. The value of a 
is arbitrary because there is no intrinsic length associated 
with the photon. 

To exhibit the geometric structure of the wave func- 
tions (55) we shall rewrite them as components of a 
Cartesian vector / = (.f xi f y ,f z ) in Cartesian coordi- 
nates, 



f(k) = Aa 2 (aky- 1 exp (~\{ak) 2 ) 
k x (to x k) + ik(n x k) 
In x fel 



(57) 



The presence of the unit vector in the direction n x k 
means that there is a vortex line in momentum space 
along the n direction with unit intensity. To obtain the 
formulas (55) we must choose the direction of n as the z 
axis in spherical coordinates. 

The increase of the lower bound in the uncertainty re- 
lation (3) from the value 3ft/2 underscores the unique 
properties of photons. This increase is due to the spe- 
cific angular dependence of the photon wave function in 
momentum space enforced by the nontrivial geometry on 
the light cone. As a result, all three functions (55) vanish 
at k = 0, in contrast to the Gaussian functions saturat- 
ing the standard Heisenberg relation. This effect was 
also present in our previous photon uncertainty relation 
[3]. In both cases, the angular dependence is the same. 
However, the radial dependence is different and this dif- 
ference is reflected in the values of the lower bounds. The 
scaling of n is chosen so that (as in [3]) the uncertainties 
in position and momentum are equally distributed, 

(R- R) = a 2 {I + \/5/2), (58a) 
(P-P) = (h/a) 2 (l + V5/2). (58b) 

Of course, their product is scale independent and gives 
the lower bound. 



VI. UNCERTAINTY RELATION FOR THE 
PRODUCT OF AR 2 AND AP 2 

The information gained in the analysis of the uncer- 
tainty relation (3) will now be used to improve the bound 
in (2). The first observation is that (R) = for all three 



'9/4+V5" 



AR 2 AP 2 








<P> 2 



FIG. 1. Dependence of the product of variances (in units 
of h 2 ) on the squared mean momentum (in units of h 2 /a 2 ). 
The leftmost dot represents the exact value of j 2 obtained for 
vanishing average momentum. The remaining dots (from left 
to right) mark the values obtained for the trial functions (60) 
with none, one, two, and up to six terms. The solid curve 
represents the two-parameter fit (61). 



functions that saturate (3), whereas the value of (P) does 
not vanish for the states with m = ±1, 



(P) 



rt3/2 + V5/2)ft +aaza K 

± -zr 1 — -- n = ±0.686- n. 

2r(l + V5/2) a a 



(59) 



Thus, already in this simple case the value of AR 2 AP 2 = 
(R ■ R)(P ■ P) - (R- R)(P) 2 is lower than the value of 
(R-R){P-P). 

In the general case, the bigger (P) is, the bigger will 
be the average photon energy. Thus, the noncommu- 
tativity of the components of R plays a decreasing role, 
bringing us closer to the situation in nonrelativistic quan- 
tum mechanics. This is clearly seen in Fig. 1 where we 
show the exact value (3) obtained for the vanishing mean 
momentum and the results of numerical calculations of 
AR 2 AP 2 . The points in this plot were obtained by 
choosing the trial functions as the product of f± (fc, 9, (f>) 
and a polynomial in k cos 9, 



1 + aikcos9 + a 2 (fccos#) 2 + a 3 (k cos 9) 3 



(60) 



where are variational parameters. These parameters 
are determined by requiring that they give the lowest 
value of AR 2 AP 2 . The points in Fig. 1 represent the 
values obtained with none, one, two, and up to six param- 
eters. The solid line represents a simple two-parameter 
fit of the form 



9/4 + 



V5 



1 + 1.14 (P) 2 + 0.8 (P) 4 ' 



(61) 



to all eight results. The numerical results clearly show 
the convergence to the value 9/4 when (P) tends to in- 
finity. This result is also in agreement with the formula 
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(13) for the commutator of the center-of-energy operators 
since the right-hand side tends to zero with the increase 
of the energy, so that at infinite energy these operators 
behave as their nonrelativistic counterparts. We shall 
confirm now this result with analytic considerations. We 
show that in the infinite-momentum frame, we indeed 
obtain as a lower bound in the uncertainty relation the 
limiting value 37i/2. Thus, our aim is to find the minimal 
value of the expression 



AR 2 AP 2 = ({A- (-R)) 2 ) ((p-(P) 



(62) 



in the limit of infinite (P) . 

In the first step, we eliminate (R) by applying the uni- 
tary transformation exp(— i(R) ■ P/h) (i.e., by choosing 
the center of the energy as the origin of the coordinate 
system). The elimination of (P) by the unitary trans- 
formation exp(i(P) • R/ft) is not so painless because the 
components of R do not commute and we are left with 
the expression 



1 



2 _ J_ l e -i(P)-R/h R.R e i(P)-R/H 



P P 



(63) 



which is to be minimized [12]. In Appendix B we find 
by the variational procedure that the minimum of 7 is 
indeed equal to 3H/2. 

The significant simplification of relativistic dynamics 
in the infinite momentum frame was noted a long time 
ago [14]. In particular, it has been shown [15] that in this 
limit, the symmetry group in the transverse plane is the 
Galilean group in two dimensions that governs nonrel- 
ativistic quantum mechanics. This explains why a(fe), 
given in the infinite- momentum frame by (Bl), is very 
simple leading to the nonrelativistic lower bound in the 
uncertainty relation for photons (2). 



VII. UNCERTAINTY RELATION FOR 
PHOTON BEAMS 

In most experiments photons appear in the form of 
photon beams. In this section we derive the uncertainty 
relation for a very common representation of such beams: 



the coherent state of the electromagnetic field. The exact 
determination of the uncertainty relation for the coher- 
ent state does not seem to be feasible but the important 
case — the limit when the mean photon number (N) is 
large — is tractable. 

Coherent states |coh) are generated from the vacuum 
state by the unitary Glauber displacement operator D 
[16], 



VW)%\ [h(k)a\{k) f* x (k)a x (k) 



D = exp 



|coh) = D|0), (64) 

where the function f\(k) that so far represented a single- 
photon state now describes an arbitrary nonmonochro- 
matic mode of electromagnetic radiation [17]. We pulled 
out the square root of the mean photon number (N) in 
the coherent state to have better control of the large (N) 
limit. The function / will be normalized to 1 as in (27). 

Our aim, as in Sec. VI, is to minimize the left-hand 
side of the uncertainty relation (62). This time [12], all 
expectation values arc to be evaluated in the coherent 
state (64), 
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((R-R)-(R)-(R)) ((P.P)-(P).(P)) 



h 2 



(65) 



A fairly complicated evaluation of the two factors appear- 
ing in this formula is relegated to Appendix C. Using the 
formulas (C8) and (C6) we obtain the following expres- 
sion for 7 2 valid for large values of (N) : 



7 = 



^k\D x f x {k)\ 2 ^k\f x {kf 



(66) 



In what follows we will tacitly assume that all results are 
valid only in the limit when (N) — > 00, and we will omit 
the symbol 0{1/{N)). 

Before subjecting this expression to the variational 
procedure, let us note that it does not depend on the nor- 
malization of f X (k). Therefore, we may vary the function 
f x (k) freely, as we did in all previous cases. The varia- 
tion with respect to f x (k) leads to the following equation 
for f x in the spherical coordinate system: 



— d K K^d K — ; 

k sin 1 



(dg sin Odg + d 2 



A 2 2-i 2 
A 2 - 2i\cos9 du) - — 

K K 



(67) 



We omitted here all intermediate steps because they are parameter k is defined now as 
analogous to those followed in Sec. V. The dimensionless 

_ hk (H) 



t\h(k)\ 2 

— . (68) 

^k\f x (k)\ 2 
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After the separation of variable we obtain the following 
equation for the radial part: 



-—d K K 3 d K + 



J(J + 1) - A 2 



K 



K{k) = — 7 2 /C(k), 
(69) 



while the angular part is the same as in (49), so that the 
lowest allowed value of j is 1. The equation for the radial 
part after the substitution, 



JC{k) = «yj'(J+i)-i cxp (_ 7 «) k,{k), 



(70) 



reduces to the equation for the confluent hypergeometric 
function, 

£(«)= iFi(|-7, W 2 7 k), (71) 



where /i = 1 + + 1)- To obtain a regular solu- 

tion, i_Fi must become a polynomial and this leads to 
the quantization condition for the parameter 7 , 



T _2»±i±£2S±iI, „ = ,1,2.... ,72, 
The lowest value of 7 is obtained for j = 1 and n = 0, 
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(73) 



Again, as in Sec. VI, the solution corresponding to the 
lowest value of 7 has a threefold degeneracy. The three 
normalized solutions, which are the counterparts of (55), 
are 



fo{k,6,<t>) = Aasm8 (ak) 



V2-l-~{ak 



f±(k,6,<f>)=Aa 



where 



(1 ± cos 6) p ± i<t> ^ ak ^-i p - iak 



V2 



(27) 



V2 



8tiT(2\/2) ' 



(74a) 
\ (74b) 



(75) 



and the parameter a sets the scale as in the case of a 
single photon. 



VIII. OBSERVABLE CONSEQUENCES OF 
UNCERTAINTY RELATIONS 

Our uncertainty relations for individual photons can be 
connected with observations through the Glauber theory 
of photodetection [18], as we have indicated in [3]. The 
interpretation of the uncertainty relation for photons is 
basically the same as in the case of the standard Heisen- 
berg uncertainty relation. The only difference is that the 
photodetection relies on the energy density of photons — 
the photon is where its energy is localized — rather than 



on the probability density to find the particle (its charge 
or mass) at a given location. To test our uncertainty re- 
lation, one would have to make repeated measurements 
on photons produced by the same source. 

The uncertainty relation plays a different role in the 
case of photon beams. In this case, the limitation on 
the dispersion AR 2 imposed by the uncertainty relation 
finds its physical interpretation in terms of the directly 
observable quantity: the focal volume. Of course, the 
focal volume does not have sharp boundaries. However, 
the moments of the energy distribution give reasonable 
measures of its size. Thus, a sensible measure of the size 
of the focal volume Vf is: 



V f = (AR 2 f\ 



(76) 



The uncertainty relation in three dimensions gives precise 
bounds on the size of the focal volume for a given spectral 
composition of the beam. According to this relation, the 
decrease of Vf is limited by the dispersion of momentum: 



V f > 



ft 3 7 3 



(AP 2 ) 



3/2 • 



(77) 



It is worth mentioning here that the one-dimensional un- 
certainty relation (15) can give only a rough estimate of 
the focal volume due to the strong correlations imposed 
by the noncommutativity of the components of R. 

The reduction of the size of the focal volume is im- 
portant in many practical applications of laser beams, 
such as fluorescence microscopy, optical tweezers, mate- 
rial processing and also in medicine. We are far from 
suggesting that our uncertainty relations will lead to an 
improvement in any of these techniques, but we believe 
that they are relevant at the fundamental level. 



IX. CONCLUSIONS 

In this work, we based the uncertainty relation for pho- 
tons on a measure of the spatial extension of the photon 
wave function, which is built around the center-of-energy 
vector: the first moment of the energy distribution di- 
vided by the total energy. By replacing the second mo- 
ment of energy used in Ref. [3] by the first moment of 
energy, we were able to bring the analysis closer to the 
standard quantum-mechanical treatment. 

The center-of-energy vector turned out to be a very 
good substitute for the nonexistent photon position op- 
erator, although the noncommutativity of its components 
leads to significant differences compared to the nonrela- 
tivistic case. In nonrelativisti c Hei s enberg uncertainty 
relations, the lowest value of V Ail 2 V AP 2 does not de- 
pend on the average position and on the average momen- 
tum. It_is not so for photons. The lowest possible value 
of V AP 2 \/AP 2 depends on the choice of the Lorentz 
frame. It varies between 3/2/1(1 + 4y/E/9) and 3ft/2, 
when the average momentum changes from to infinity. 
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Somewhat paradoxically, highly energetic photons obey 
almost the same uncertainty relations as nonrelativistic 
particles. This is explained by the special properties of 
rclativistic dynamics in the infinite-momentum frame. 

The uncertainty relations based on the center-of- 
energy operator were also derived for photon beams de- 
scribed by coherent states of the electromagnetic field. 
Analytic results were obtained in the limit of a large num- 
ber of photons in the beam. These uncertainty relation 
give a fundamental limitation on the reduction of the 
beam focal volume. 



Finally, using the lemma, we may replace H by y H in 
the first term and expand the resulting double commu- 
tator: 



= 



= —N - 
H 



N 



H 



Nl 
H 



N 



H 



H 



(A5) 
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The vanishing of the difference of two expressions for R 
appearing in (9) means that they are equal. 



Appendix A 

To prove the equality of the two forms of R in (9) we 
will first prove the following lemma: 



If [H, C]=0 then [VH, C] = 0. 



(Al) 



In the proof, we use the fact that the eigenvectors of the 
Hamiltonian form a basis. Acting on an arbitrary state 
in this basis \E) (excluding the vacuum), we have 

(yi+y/Ey^/I,C]\E) = [H,C]\E)=0. (A2) 

Since the factor does not vanish, it can be 

dropped and the validity of the lemma is established. 

Next, we use the commutation relations between the 
Hamiltonian and the generator of the Lorentz transfor- 
mations 



[H,N] 



-ihP, 



(A3) 



Appendix B 



To apply the variational procedure, we rewrite the 
functional (63) in the one-photon space. To simplify the 
calculations, we choose n in the direction of the aver- 
age momentum. With this choice, the operator in-R re- 
duces to an ordinary derivative with respect to k z because 
the scalar product n ot vanishes. Therefore, the unitary 
operator e l ( p )- R / h acting on the photon wave functions 
becomes just the shift operator. Therefore, in the func- 
tional (39) the argument k z of a(k) is shifted now by 
(k z ) = (P z /h). The solution of the differential equation 
obtained by varying is a very difficult task because 
the variables k and 9 can no longer be separated. 

However, in the limiting case when (k z ) tends to in- 
finity, there is a radical simplification. In this limit, 
a(k + (fe)) becomes 



to obtain 



1 N 1 



H 

H,N 



H 
1 



H V H 



P 



o. 

(A4) 



D y 



and the variational equation (45) is replaced now by 
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1 



k 2 sin 2 9 



K - 2 7 



which again allows for the separation of variables, 



g x {K,6,<p)=K{K)Q(6)e im v. 



(B3) 



The radial and the angular parts satisfy the equations 



K K 



■ ae sin Odg + 



K(n) = 2 7 /C(k), (B4) 



sin( 



sin 2 6 



0(6) =j(j + 1)0(6). 

(B5) 



This time, the value j = is allowed provided we choose 
m = A. The equation for the radial part is that of the 
spherically symmetric harmonic oscillator. The lowest 
value of 7 = 3/1/2 is obtained for the ground state. This 

I 



g\(K,6,ip) = 0, 



(B2) 



r 



confirms the prediction made on the basis of our numer- 
ical calculations. 



Appendix C 

In this appendix, we evaluate the leading terms of the 
expansion in 1/(N) for the dispersion in position (21) 
and momentum (22) evaluated in the coherent state of 
the electromagnetic field. In this calculation, we use the 
second form of the position operator (9). The expectation 
value of any combination of creation and annihilation op- 
erators in a coherent state is tantamount to the vacuum 
expectation value of the same combination of these op- 
erators transformed by the action of the displacement 
operator D, 

rta[(k)D = a[(k) + ^{N)r x (k), (Cla) 
D^ x (k)D = ax(k) + y/{Njfx{k). (Clb) 

We will need only the following lowest-order correction 
to the operators H, N, and P: 



D^HD = (N}hc 
D^ND = (N)hc 
D^PD = (N)h 



^ft(k)fx(k) + -^=^(a{(k)f x (k) + f* x (k)a x (k)) + 0(^L) 
^fl(k)iD x f x (k) + -^^(a{(k)iD x f x (k) + r x (k)iDxa x (k)) +0(j^) 
^f{(k)nf x (n) + -^^(a[(k)nf x (k) + f* x (k)na x (k)) +0(j^) 



The first two formulas lead to the following expression for R: 

1 " 



n 



W) 



^(a[(k)f x (k) + r x (k)a x (k)) +0(j^) 



(C2a) 
(C2b) 
(C2c) 



(C3) 



where 



u 



^f* x (k)h(k), (C4a) 

Ar = ^f* x (k)iDxfx(k) (C4b) 
K = Af/H. (C4c) 



I 

In both factors of (65), the leading terms cancel because 
they are c numbers, so that there is no difference be- 
tween the averaged square and the square of the average. 
We shall first calculate the next-order corrections to the 
difference (R-R) — (R)-(R). First, note that if the con- 
tribution comes from only one R, then it does not con- 
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tribute to the difference because it cancels out between 
the two terms. The 0(1/ (N)) terms are not canceled by 
their counterparts in (R) ■ (R) only when the corrections 
appear in both operators R in R R. The same obser- 
vation holds for the momentum operator. Therefore, the 
lowest-order corrections come only from the products of 
two terms linear in the creation and annihilation opera- 



tors, and the final results can be written in the form 

<(*-^) 2 )=^W 

x ^ k [{iD x - H)f x (k)]*-(iD x - n)h(k), (C5) 

^(p- (P>)^ = (N)h^kf*(k)f x (k). (C6) 

Without any loss of generality [the function f\(k) is at 
this point arbitrary and it will be determined from the 
variational procedure later], we can make the following 
replacement: 

h(k) -+ exp(-ife-^)/ A (fc). (C7) 

This change of phase makes no difference in (C6), but 
it leads to the elimination of the 7?.- dependent terms in 
(C5), and we obtain 

((A-<fl>) a ) = ^y^*|I?A/A(*)| a . (C8) 
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